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Feynman's "no-node" theorem states that the conventional many-body ground-state 
wavefunctions of bosons in the coordinate representation is positive-definite. This im- 
plies that time-reversal symmetry cannot be spontaneously broken. In this article, we 
review our progress in studying a class of new states of unconventional Bose-Einstein 
condensations beyond this paradigm. These states can either be the long-lived metal- 
stable states of ultra-cold bosons in high orbital bands in optical lattices as a result of 
, the "orbital-Hund's rule" interaction, or the ground states of spinful bosons with spin- 

^ ' orbit coupling linearly dependent on momentum. In both cases, Feynman's argument 

' does not apply. The resultant many-body wavefunctions are complex-valued and thus 

break time-reversal symmetry spontaneously. Exotic phenomena in these states include 
the Bose-Einstein condensation at non-zero momentum, the ordering of orbital angular 
momentum moments, the half-quantum vortex, and the spin texture of skyrmions. 

Keywords: Bose-Einstein condensation, optical lattices, exciton, time-reversal symmetry, 
^ ~ ' spontaneous symmetry breaking 
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5-H ' ground state wavefunctions of bosons are positive- definite in the coordinate repre- 

sentation provided no external rotation is applied and interactions are short-ranged 
The proof is very intuitive: due to the time-reversal (TR) symmetry, the ground 
state wavefunction Vl/(ri, ....,fvi) can be chosen as real. If it is not positive-definite, 
i.e., it has nodes, the following surgery can be done to lower its energy. We first 
take its absolute value of |^'(ri, r„)|, whose energy expectation value is exactly 
the same as that of '^{ri, ....,r„). However, such a wavefunction has kinks at node 
points. Further smoothing the kinks results in a positive-definite wavefunction, 
and the kinetic energy is lowered by softening the gradients of the wavefunction. 
Although the single body potential energy and the two-body interaction energy 
increase, they are small costs of a high order compared to the gain of the kinetic 
energy. We can further conclude that the ground state is non-degenerate because 



1. Introduction 

In Feynman's statistical mechanics textbook, it is stated that the many-body 
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two degenerate positive-definite wavefunctions cannot be orthogonal to each other. 

This "no-node" theory is a very general statement, which applies to almost all 
the well-known ground states of bosons, including the superfluid, Mott-insulating, 
density-wave, and even super-solid ground states. It is also a very strong statement, 
which reduces the generally speaking complex-valued many-body wavefunctions to 
be positive-definite distributions. This is why the ground state properties of bosonic 
systems, such as "^He, can in principle be exactly simulated by the quantum Monte- 
Carlo method free of the sign problem. Furthermore, this statement implies that the 
ordinary ground states of bosons, including Bose-Einstein condensations (BEC) and 
Mott-insulating states, cannot spontaneously break time-reversal (TR) symmetry, 
since TR transformation for the single component bosons is simply the operation 
of the complex conjugation. 

It would be exciting to search for exotic emergent states of bosons beyond this 
"no-node" paradigm, whose many-body wavefunctions can be complex-valued with 
spontaneous TR symmetry breaking. Since properties of complex-valued functions 
are much richer than those of real- valued ones, we expect such states can exhibit 
more intriguing properties than the ordinary ground states of bosons. For this 
purpose, we have recently made much progress with two different ways to bypass 
Feynman's argument, i ncludin g the meta-stable states of bosons in the high orbital 
bands of optical lattices ^^EIll and multi-component bosons with spin-orbit coupling 
linearly dependent on momentum ^ . 

Clearly, the "no-node" theorem is a ground state property which does not apply 
to the excited states of bosons. The recent rapid development of optical lattices with 
ultra-cold bosons provides a wonderful opportunity to investigate the meta-stable 
states of bosons pumped into high orbital bands. Due to the lack of dissipation 
channels, the life time can be long enough to develop inter-site coherence " ' . We 
have shown that the interaction among orbital bosons are characterized by the 
"orbital Hund's rule" which gives rise to a class of complex superfluid states by 
developing the on-site orbital angular momentum (0AM) moments. In the lattice 
systems, the inter-site hoppings of bosons lock the 0AM polarization into regular 
patterns depending the concrete lattice structures. 

Furthermore, the "no-node" theorem even does not apply to the ground states 
of bosons if their Hamiltonians linearly depend on momentum, i.e., the gradient op- 
erator. A trivial example is the formation of vortices in Bose condensates with the 
external rotation in which the Coriolis force is represented as the vector potential 
linearly coupled to momentum. However, in this case TR symmetry is explicitly 
broken. A non-trivial example is that bosons with spin-orbit (SO) coupling, whose 
Hamiltonians also linearly depend on momentum and are TR invariant. The inva- 
lidity of the "no-node" theorem can also give rise to complex-valued ground state 
wavefunctions. 

Although "^He is spinless and most bosonic alkali atoms are too heavy to exhibit 
the relativistic SO coupling in their center-of-mass motion, SO coupling can be 
important in exciton systems in two dimensional quantum wells. We have shown 
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that the Rashba SO couphng in the conduction electron bands induces the same 
type of SO coupling in the center-of-mass motion of excitons. In a harmonic trap, 
the ground state condensate wavefunction can spontaneously develop half-quantum 
vortex structure and the skyrmion type of the spin texture configuration'^. On the 
other hand, effective SO coupling in boson systems can be induced by las er beams, 
which has been investigated in several publications by other groups in 

RefsmHH 

In the following, we will review our work in both directions outlined above 
including many new results unpublished before. In Section [21 we explain the char- 
acteristic feature of interacting bosons in high orbital bands, the "orbital Hund's 
rule" , and the consequential complex-superfluid states with the ordering of on-site 
0AM moments. The ordering of 0AM moments in the Mott-insulating states is 
also investigated. In Section[3l we review the TR symmetry breaking states of spin- 
ful bosons with spin-orbit coupling. Interesting properties including half-quantum 
vortex, and the skyrmion-like spin-textures are studied. Conclusions are made in 
Section m 

Due to the limit of space, we will not cover interesting related topics of orbital 
bosons, such as the nematic superfluid state ^andthe algebraic superfluid state^^, 
and the orbital physics with cold fermions | l'^ | 14 | 15[r6 | 17 | 18 | 19 ( 20 |^ which has also 
aroused much research attention recently. 



2. "Complex-condensation" of bosons in high orbital bands in 
optical lattices 

In this section, we will review the "complex-condensation" with TR symmetry 
breaking, which is a new state of the p-orbital bosons. This is an example of novel 
orbital physics in optical lattice with cold bosons. Below let us give a brief general 
introduction to orbital physics for the general audience. 

Orbital is a degree of freedom independent of charge and spin. It plays im- 
port antroleshi magnetism, superconductivity, and transport in transition metal 
oxides 12222123] 'yy^q key features of orbital physics are orbital degeneracy and spa- 
tial anisotropy. Optical lattices bring new features to orbital physics which are not 
easily accessible in solid state orbital systems. First, optical lattices are rigid lat- 
tices and free from the Jahn- Teller distortion, thus orbital degeneracy is robust. 
Second, the meta-stable bosons pumped into hi gh orbital bands exhibit novel su- 
perfluidity beyond Feynman's "no-node" theory^ ^ l ^^ l '^ l ^^ l ^^ l '^ i Third, p-orbitals 
have stronger spatial anisotropy than that of d and /-orbitals, while correlation 
effects in p-orbital solid state systems (e.g. semiconductors) are not that strong. In 
contrast, interaction strength in optical lattices is tunable. We can integrate strong 

correlation with spatial anisotropy more closely than ever in p-orbital optical lattice 
systems ^3 14 15 16 17_ 

Rece ntly, orbital ^ has been attracting a great deal of 

attentioii^ ^ l ^ l ^^ l ^ l ^^ l ^ ^'l^'^'"^^"^^"^'^'^-^. For orbital bosons, a series of theory 
works have been done [6 | 26 | 2 | 24 | 3 , 12 , 25 i 4l^ ^ including the illustration of the ferro- 
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Fig. 1. A single site problem with two spinless bosons occupying p^j^y-orbitals: (A) the OAM singlet 
and (B) one of the OAM doublets. The latter is energetically more favorable than the former as 
a result of the "orbital Hund's rule" . The bonding pattern of p-orbitals: (C) the tr-bonding and 
(D) the TT-bonding. 



orbital nature of interac tion s^, orbital superfluidity with spontaneous time-reversal 
symmetry breaking 24l4l^ ^^^^ nematic superfluidity- xhe theory work on the 
p-orbital fermions is also exciting includi ng the fla. t band and associated strong cor- 
relation physics in the honeycom b lattic el '*^'^ ! ^^ ! ^ -^, orbital exchange and frustration 
| 16 | 20| g^^^ topological insulators tSGH 

The experimental progress has been truly exciting |29 | 30 | 31 [7l Mueller et al. 
^ have realized the meta-stable p-orbital boson systems by using the stimulated 
Raman transition. The spatially anisotropic phase coherence patt ern has been ob- 
served in the time of flight experiments. Sebby-strabley et al. have successfully 
pumped bosons into excited bands in the double-well lattice. In addition, p-orbital 
Bose-Einstein condensation (BE C) h as also been observed in quasi one-dimensional 
exciton-polariton lattice systems 

Below we illustrate the important feature of interactions between orbital bosons, 
the "orbital Hund's rule" , which results in TR symmetry breaking. 



2.1. Orbital Hund's rule of interacting orbital bosons 

The most remarkable feature of interacting bosons in high orbital bands is that 
they f avor to maximize their onsite orbital angular momentum (OAM) as explained 

below 1213!. 

Let us illustrate this Hund's rule type physics through the simplest example: a 
single site problem with two degenerate p^^^j^-orbitals filled with two spinless bosons. 
Because bosons are indistinguishable, the Hilbert space for the two-body states only 
contains three states, which can be classified according to their OAM: an OAM 
singlet as ^{plpl + pIpI)\0) depicted in Fig [T] A, and a pair of OAM doublets 
■^^{pI ± *Py)^|0) with Lz — ±2 as depicted in Fig[T]B. Assuming a contact inter- 
action V = gS(ri — r2), the interaction energy of the former is calculated as 
while that of the latter is |[/ with the definition of 
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Fig. 2. The p-orbital band structure in the two dimensions: A) square and B) triangular lattices. 
In the square lattice, the band minima are located at Kx = {7r/ao, 0) for the pa;-orbital band, and 
Ky = (0, 7r/ao) for the py-orbital band, respectively. In the triangular lattice, the band minima 
are at three middle points of Brillouin zone edges as Ki = (0, -^y-) and i4'2,3 = (-^,111 ^ ). The 
corresponding orbital configurations at these three band minimal are polar-like and parallel to the 
momentum directions of Ki^2,3 respectively as depicted. 



In the 0AM singlet state, bosons occupy polar (real) orbitals (e.g. px,Py), whose 
angular distribution in real space is narrower than that of the axial (complex) 
orbitals (e.g. Px ± ipy ) for the 0AM doublets. By occupying the same axial orbital 
and therefore maximizing 0AM, two bosons enjoy more room to avoid each other. 

This "ferro-orbital" interaction is captured by the following multi-band Hubbard 
Hamiltonian for the p-orbital bosons as 



where n = p^Px + p\Py and = —i[p\py — p\px)- The first term is the ordi- 
nary Hubbard interaction and the second term arises because of the orbital degree 
of freedom. In three dimensional systems in which all of the three p-orbitals are 
present, we only need to replace with = + Lif^ + L\ and Lx^y defined as 
i^x = -i{v\Vz - pIpv), Ly = -i{p\px - pIpz)- 

When more than two bosons occupy a single site, bosons prefer to go to the 
same single particle state by their statistical properties. Again going into the same 
axial state and thus maximizing 0AM can reduce their repulsive interaction energy. 
This is an analogy to the Hund's rule of electron's filling in atomic shells. The 
first Hund's rule of electrons maximizes electron total spin to antisymmetrize their 
wavefunction, and the second Hund's rule further maximize their 0AM to extend 
the spatial volume of wavefunction. The key feature is that electrons want to avoid 
each other as possible as they can. For the orbital physics of spinless bosons, the 
same spirit applies with the maximizing of 0AM. 

Now we have understood the single site physics in which bosons develop rotation. 
From the symmetry point of view, it is similar to the Px + ipy superconductors. 
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2.2. Band structures of the p-orbital systems 

The p-orbital Hamiltonian within the tight-binding approximation can be written 

as 

where the unit vector Cij is along the bond orientation between two neighboring sites 
i and j and /y = z x Bij is perpendicular to e^j. Pe^j and p^ are the projections of 
p-orbitals along (perpendicular to) the bond direction respectively as defined below 

= {Pxex +Pyey) ■ Pf^^ = {PxS-x + Py^y) ' fir (4) 

The cr-bonding t|| and the 7r-bonding t± describe the hoppings along and perpendic- 
ular to the bond direction as depicted in Fig.[T]C and D, respectively. The opposite 
signs of the a and 7r-bondings are due to the odd parity of p-orbitals. t± is usually 
much smaller than in strong periodical potentials. 

E g. B ] exhibits degenerate band minima in both the square and triangular 
lattictpEGSl. In the square lattice, the Brillouin zone (BZ) is a square with the 
edge length of where qq is the length between the nearest neighbors. The spec- 
tra read ep^ = i|| cosfcj; — t^cosky and = —t±coskx + cos ky. As depicted 
in Fig. [21 A, the band minima are located at Kp^ — (^,0) for the pa;-band and 
Kpy = (0, — ) for the Pj,-band, respectively. 

In the triangular lattice, we set the 7r-bonding to zero which does not change the 
qualitative band structures. We take unit vectors from one site to its six neighbors 
as ±ei.2.3 with ii — Bx, B23 — —\Bx ± The Brillouin zone takes the shape of 

a regular hexagon with the edge length 47r/(3a). The energy spectrum of Hq is 

^W=i||{.4T//|-35,r}, (5) 

with 

3 

/g = ^cos(fc • Ci), 5fc= X! cos(fc • ei)cos(^- Cj). (6) 

i=l 3>j:>j>l 

The spectrum contains three degenerate minima located at the three non-equivalent 
middle points of the edges as = (0, 7^), ^2,3 = (±f , The factor e^^i *" 

takes the value of ±1 uniformly in each horizontal row but alternating in adjacent 
rows. If the above pattern is rotated at angles of then we arrive at the patterns 
of e*^^'^ ''. Each eigenvector is a 2-component superposition vector of px and Py 
orbitals. The eigenvectors at energy minima are ipKi = ^^^^'^\Py)- i^K2.3 can be 
obtained by rotating -0i at angles of ±^ respectively as depicted in Fig. [2]B. 

2.3. Weak coupling analysis - BEC at nonzero momenta 

If there were no interactions, any linear superposition of the above band minima 
would be a valid condensate wavefunctions. However, interactions will select a par- 
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B) 



Px+'tV 




(units of n/a) 



Px-ip 



Fig. 3. A) The time-of-flight spectra of the p-orbital Bose condensation in the square lattice. The 
coherence peaks occurs at (m^, (n + ^)^), and ((m + B) The staggered ordering 

of OAM moments in the square lattice with the phase pattern on each site. Each of the a and 
TT-bonds achieves phase coherence. From Liu and Wu, Ref 2. 

ticular type combination which exhibits TR symmetry breaking and agrees with 
the above picture of "orbital Hund's rule" . 

2.3.1. The square lattice - the staggered ordering of OAM moments 

In the square lattice, we take the condensate wavefunction as ^sgC^O = 
ci'4}p^{Kp^) + C2ippy{Kp^) under the constraint of |cip + |c2p = 1. Any choice 
of Ci^2 minimizes the kinetic energy. However, the interaction [/-terms break the 
degeneracy. The degenerate perturbation theory shows that the ground state values 
of Ci^2 take Ci = 1, C2 = « or its equivalent TR partner of Ci — 1, C2 = —i. The mean 
field condensate can be described as 



where Nq is the particle number in the condensate. 

For better insight, we transform the above momentum space condensate to the 
real space. The orbital configuration on each site reads 



where the U{1) phase e"^** is specified at the right lobe of the p-orbital. The Ising 
variable = ±1 denotes the direction of the OAM, and is represented by the anti- 
clockwise (clockwise) arrow on each site in Fig. |3]B. Each site exhibits a nonzero 
OAM moment and breaks TR symmetry. The condensate wavefunction of Eq. [7] 
describes the staggered ordering of af. We check that the phase difference is zero 
along each bond, and thus no inter-site bond current exists. 

The condensate described in Eq. [7] breaks both TR reversal and translational 
symmetries, and thus corresponds to a BEG at non-zero momenta. This feature 
should exhibit itself in the time-of-flight experiments which have been widely used 
to probe the momentum distribution of cold atoms. The coherence peaks of Eq.[7]are 




e^'^ilPx) +i<Jf\Pv)), 



(8) 
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not located at integer values of the reciprocal lattice vectors but at (^ti^, ("•+5)^)) 
and ((m + as depicted in Fig. [3] A. Furthermore unlike the s-orbital 

condensate the p-wave Wannier function superposes a non-trivial profile on the 
height of density peaks. As a result, the highest peaks are shifted from the origin — 
a standard for the s-wave peak — to the reciprocal lattice vectors whose magnitude 
is around \/l where / is the characteristic length scale of the harmonic potential of 
each optical site. A detailed calculation of form factors is given in Ref. 



2.3.2. Excitations - gapless phonons and gapped orbital modes 

The elementary excitations in p-orbital BECs consist of both the gapless phonon 
mode and the gapped orbital mode. The former corresponds to the Goldstone mode 
related to the U{1) symmetry breaking, and the latter corresponds to the flipping 
of the direction of 0AM moments. In the following, we will take the 2D square 
lattice as an example. 

We assume the condensate as -^[pHQx) +'ipl{Qy)] with ~ (tt, 0) and Qy = 
(0,7r). The boson operators take the expectation values as 

(b,(f)|> = {-ly^cb., {\py{r)\) ^ z(-ir«<^. (9) 

From minimizing the onsite part of the free energy respect to the condensate order 
parameter (j) 

F = -,,n+^{n'-hl), (10) 

we have /i = respect to the band minima. The fluctuation around the ex- 

pectation value is defined as 

Px{r) = {\Pxir)) + Spx, Py{r) = {\py{r)) + 5py. (11) 

Then the interaction terms in Eq. [2] are expanded as 
4 8 1 

Hint = g 101'* + [SPISPX + 5pl5py\-\(j)\^ + -(I)*''^[5px5px - SpySpy] 

+ lc^'[6plSpi - SpISpI] + ?^{-y^+ry[cj,HplSpl - cj^*^'6px5py]. (12) 
Combined with the free part, we arrive the mean field Hamiltonian as 

HmF = J2 *^(^)aM„b(fc)*,(fc) (13) 

k,a,b 

where '^/'^{k) = [pUk),pl{k + Q),p^{-k),py{~k - Q)], Q = (7r,7r) = - Qy, and 
the summation is only over half of the Brillouin zone. The matrix kernel reads 

/e.(^) + f|0p f</.2 ,^^2 ^^ 

e^(fc + Q) + 



M(fc) 



V -fc/r-' fc/r'' e^(fc + g) + i|i|^|2y 
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tried 






strong couplingweak coupling 
a=7i/4 a=it/6 



-2-1 1 2 3 
kj (units of 2n/a) 



Fig. 4. A) The time-of-fiight spectra of the p-orbital boson condensation in the triangular lattice. 
The coherence peaks occurs at B) The stripe ordering of OAM moments in the triangular lattice. 
From Wu et al Ref. 3. 



The spectra are reduced to generalized eigenvalue problem of 

X-i{diag(l, 1, -1, -1) M{k)]X = diag(Si, ^2, -Ei), 



(14) 



where i5i,2,3,4(fc) are excitation eigenvalues, and X contains the eigenvectors. 

Let us consider the excitation close to the condensation wavevector k = Qx + q, 
and then k + Q ~ Qy + q. At small value of q, we obtain the excitation spectra as 



EiAk)~- 

wheree(g) = ^[ex{Qx + q) + £y{Qy + q)-ex{Qx)-ey{Qy)] ~' 



(15) 

-{ql + qD- Clearly, 



the gapless mode with the linear dispersion relation describes the superfluid phase 
fluctuations; the gapped mode describes the orbital excitations corresponding to 
the flipping of orbital angular momenta. 



2.3.3. The triangular lattice - the stripe ordering of OAM moments 

In the triangular lattice, the OAM moments instead form a stripe ordering, i.e., 
the OAM moments along one row polarizes along the z-axis and those along the 
neighouring rows polarizes with the opposite direction. This can be intuitively un- 
derstood as follows. In the superfluid state, the OAM moments behave like vortices 
whose interactions are long range. The above stripe configuration of positive and 
negative vortices is the optimal configuration to minimize the globe vorticity. 

Let us first examine the weak coupling limit. Again we write a general form for 
condensation wavefunction as a linear superposition of the three band minima 

= -y={ci'ipKi + C2ipK2 + CsTpKa), (16) 

V N 

where -^^1,2,3 are the locations of band minima defined in Sect. 12.21 Without loss of 
generality, we set ci = 1, and define x = (c2 + C3)/2, y = (c2 — C3)/2 and x = Xi+ix2 
and y ~ yi + ij/2, then the normalization factor N = -^/l + 2{x^ + yf + X2 + J/|). 
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The interaction energy persite is calculated as 

E^nt = ^{i^+^-yl- vlf + 82:? + Syl] + 1 

- ^{^{xiy2-X2yif + 2xl + 2yl]. (17) 

The terms in the first line are from the density-density interaction which can be 
minimized by setting a:i = yi = and X2 — ±j/2- This means only one of 02,3 
is nonzero and purely imaginary. In this case, the particle number on each site is 
uniform. The terms in the second line can also be minimized at this condition with 
a further requirement of X2 — ±2/2 = ±5- If we take C2 nonzero, then C2 = ±i. 

Thus the mean field condensate can be expressed as ^jy^i { "^(^^2 "'^ ^^ifa) | l'^) 
with |0) the vacuum state and A^o the particle number in the condensate. This state 
breaks the U{1) gauge symmetry, as well as TR and lattice rotation symmetries, 
thus the ground state manifold is U{\) ® Z2 ® ^3. This state also breaks lattice 
translation symmetry, which is, however, equivalent to suitable combinations of 
U{1) and lattice rotation operations. 

Again we transform the above momentum space condensate to the real space, 
whose orbital configuration takes 

e*"^'"'(cosa|px) + icTfjsinalpj^)) (18) 

with a = ^ as U/t 0. The general configuration of a is depicted in Fig.[4]B for 
later convenience. At U/t — > 0, Px.y are not equally populated, and the moment 
per particle is This does not fully optimize Hint which requires L^^p = i^- 

However, it fully optimizes Hq which dominates over Hint in the weak coupling 
limit. We check that the phase difference is zero along each bond, and thus no 
inter-site bond current exists. 

Interestingly, as depicted in Fig. [3]B, 0AM moments form a stripe order along 
each horizontal row. This stripe ordering in the weak coupling limit is robust at 
small values of the 7r-bonding t± because it does not change the location of the 
band minima and the corresponding eigenfunctions of "^Ki 2 3 ^t all. 

The driving force for this stripe formation in the SF regime is the kinetic energy, 
i.e., the phase coherence between bosons in each site. By contrast, the stripe forma- 
tion in high Tc cuprates is driven by the competition b etw een long range repulsion 
and the short range attraction in the interaction terms 

This stripe phase should manifest itself in the time of fiight (TOF) signal as de- 
picted in Fig.|4]A. In the superfiuid state, we assume the stripe ordering wavevector 
Ki, and the corresponding condensation wavevectors at K2,3. As a result, the TOF 
density peak position after a fight time of t is shifted from the reciprocal lattice 
vectors G as follows 

(n(r)>, (X 5] {|02(a, k)\'S\k~K2 ~ G) + |03(a, k)\'S\k-Ks - G)}, (19) 
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Fig. 5. A) Strong coupling analysis for the inter-site coupling between OAM moments on neigh- 
boring sites. The OAM moments are described by Ising variables aij. Bi and 6j are the azimuthal 
angles of the bonds relative the x-axis. B) The stripe ordering of OAM moments in the triangular 
lattice. The smallest unit is rhombic with the total vorticity of 1/3. 

where k = mf/{ht); 02,3(a,fc) is the Fourier transform of the Wannier p-orbital 
wavefunction |02,3(a)): and G = ^[m, (— m + 2n)/-\/3] with m,n integers. Thus 
Bragg peaks should occur at ^[m± i, -^(— m + 2n+ i)]. Due to the form factors 

of the p-wave Wannier orbit wavefunction \(j)2.3{a, ^)p7 the focations of the highest 
peaks is not focated at the origin but around \k\ « l/l^^y Due to the breaking of 
lattice rotation symmetry, the pattern of Bragg peaks can be rotated at angles of 

-1-221 

2.4. Strong coupling analysis - the lattice gauge theory formalism 

In this subsection, the ordering of the OAM moments in the strong coupling su- 
perfluid regime is examined. We will employ the lattice gauge theory formalism 
developed by Moore and Lee in Ref. in the context of the p + ip Josephson 
junction array systems. 

In the superfluid regime, each site i is denoted by a C/(l) variable which 
denotes the superfluid phase along the x-direction, and by an Ising variable az,i — 
±1 for the direction of the OAM moment as depicted in Fig. [5l Along a general 
direction of the azimuthal angle 9, the superfluid phase is + <^z,iO- The intersite 
a and 7r-bonding become the intersite Josephson coupling as 

H,j = -nt|i cos{4>i - - A\\{i,j)) - nt±cos{(j)., - (jy^ - A^^ii,])), (20) 

where n is the average particle number per site. The phase differences in Eq. 1201 
takes into account the geometric orientation of the bond {ij) as captured by the 
gauge fields A\\ and A±^ with the definition that A\\ = (yz,i9i — (yz,j0j, A± = az.i{Oi + 
^) — az,j{dj — f) = ^11 + ^{'^z,i + o'zj), and 9j = 9i + n arc the azimuth angles 
relative to the x-axis defined in Fig. [51 

Let us first consider the case oi t± =/= and write down an effective Hamiltonian 
for the Ising variables az ■ In order to minimize both the Josephson couplings of the 
cr and TT-bonds, we need at = —aj and thus = Aj__ij. Otherwise, if ctj ~ Uj, 
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then = ^±-ij + This costs an energy of 2t±, and leads to an effective 

antiferro-orbital interaction between the Ising variables as 



Heff = nt±_ az,iCTz,j- (21) 

ij 

Thus the antiferro-orbital ordering of 0AM moments in the square lattice remains 
valid in the strong coupling limit at non- vanishing t±^ , which enforces that both the 
a and 7r-bonds achieve phase coherence as depicted in Fig. [3] B. 

On the other hand, in the limit of a vanishing t^, the leading order effect involves 
mu ltiple site interaction around a plaquette. We follow the method described in Ref. 
^231, and perform the duality transformation for the U{1) phase variables </> under 
the background of the geometric gauge potential Ajj . We separate the contributions 
from the phonon part and the vortex part as 

m X 

+ 7rni|| ^(TO(a;)-$,)log^^^(m(x') -$',)}, (22) 

where Z^h is the phonon contribution in the vortex-free configuration; x marks the 
dual lattice site (or the plaquette index in the original lattice); m is the vortex 
charge; '^^ is the external flux through the plaquette defined as 

^--^Y.h.^3■ (23) 

id 

In the square lattice, the geometric gauge flux around each plaquette x is cal- 
culated as 

= - ^ (CTj + cr,+e, + cr»+e,. +e„ + CT-*+e J ■ (24) 

If is integer- valued, it can be absorbed by shifting the zero of the vortex charge 
as m! (x) = m(x) — $3, which remains integer-valued, thus there is no c ost for energy. 
The feature can be captured by the effective Hamiltonian in Ref. 1^3] as 

-ffe// = -A'^CTifJiCTfeCrz, (25) 

where z, j, fc, I are four sites around a plaquette centered at a;; -ft" « nt\\ is the energy 
scale of the a-bonding. This model has a sub-extensive Z^ symmetry investigated 
in Ref. i.e., flipping the sign of Oz along each row or column leave Eq. [25l 
invariant. Thus it cannot develop the ordering of 0AM at any finite temperature. 
Let us come back to the superfluid sector 0, its most relevant topological defect is 
the half-quantum vortices because can take the values of ±i. The Kosterlitz- 
Thouless transition is associated with the unbinding of half-quantum vortices. As 
a resul t, th e low temperature phase is the quasi-long range ordering of pairing of 
bosons 
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Next we move to the strong coupling theory in the triangular lattice. We start 
from the limit of tj^ — 0. The geometrical gauge flux becomes 

*x = -^(CTl+CT2 + cr3), (26) 


where 1,2,3 are three sites around a triangular plaquette. There is no way to form 
an integer flux. The smallest vorticity per plaquette is ±i which corresponds to 
either two +l's and one —1, or two — I's and one +1 to minimize the vortex core 
energy. In such a dense vortex configuration, Eq. [21] rigorously speaking does not 
apply because its validity relies on the assumption of small vortex fugacity. However, 
the structure of interactions among vortices still implies that vortices form a regular 
lattice with alternating positive and negative vorticity. The dual lattice (the center 
of the triangular plaquette) is the bipartite honeycomb lattice. It is tempting to 
assign ±i alternatively to each plaquette, but actually it is not possible due to 
the following reason. Consider a plaquette with vorticity +g, thus its three vertices 
are with two +l's and one —1. The neighboring plaquette sharing the edge with 
two +l's must have the same vorticity, and merges with the former one to form 
a rhombic plaquette with vorticity i as depicted in Fig. [5] B. Thus the ground 
state should exhibit a staggered pattern of rhombic plaquette with vorticity of ±i. 
This arrangement precisely corresponds to the stripe order of the Ising variables. 
This configuration breaks both the lattice rotational and translational symmetries, 
which is six-fold degenerate. If we turn on the small term, it results in an 
antiferromagnetic Ising coupling with nearest bonds, the stripe configuration also 
satisfies its ground state requirement. Thus we believe that t± term does not change 
the stripe configuration in the strong coupling limit either. 

Having established the stripe order for the 0AM moments, it is straightforward 
to further optimize the U{\) phase variable 0. The result is depicted in Fig. |4]B 
with (j) marked on the right lobe of the Px ± ipy orbit. Each horizontal bond has 
perfect phase match, while each of tilted bond has a phase difference of . Thus 
around each rhombic plaquette, the phase winding is x 4 = ^ and gives rise to 
the Josephson supercurrent along the bonds as 

.7 = ^sinA^, (27) 

with AO = ^ and the directions specified by arrows in Fig. [3]B. In other words, in 
addition to the stripe order of the Ising variables which corresponds to the onsite 
0AM moments, there exists a staggered plaquette bond current order. This feature 
does not appear in the square lattice The orbital plaquette current ordering has 
been studied in the stro ngly correlated electron systems, such as the circulating 
orbital current phase and the d-density wave states in the high Tc compounds 
ISni. It is amazing that in spite of very different microscopic mechanism and energy 
scales, the two completely different systems exhibit a similar phenomenology. 
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Fig. 6. A) Phase diagram based on the GMF theory in the 2x2 unit cell (see Fig.|4]B). Large scale 
Gutzwiller mean field calculations in a 30 X 30 lattice are performed to confirm the stripe ordered 
superfluid (SF) phase at points 1, 2 and 3 with (t/U, fi/U) = (0.02, 2), (0.03, 1.5) and (0.038, 2.2), 
respectively. B) The fiux $ around a rhombic plaquette v.s. t/{nU). It decays from i in the strong 
coupling limit to in the non-interaction limit. The solid line is the Gutzwiller result at n = 3, 
while the dashed line is based on the energy function Eg. I30l of the trial condensate. From Wu et. 
al Ref. 3. 



2.5. Intermediate coupling regime in the triangular lattice - 
self- consistent mean field analysis 

We have shown that in the triangular lattice the stripe ordering of the 0AM 
moments exist in both weak and strong coupling superfluid states. The orbital 
configurations are slightly different: the orbital mixing angle a defined in Eq. [18] 
equals to ^ in the weak coupling limit, while it equals to ^ in the strong coupling 
limit. This arises from the competition between the kinetic energy and the onsite 
Hund's rule. Below we will see that as goes from small to large, the stripe 

ordering remains the same with a smooth evolution of a from ^ to ^ as depicted 
in Fig. HE. 

We use a Gutzwiller type mean field theory for a 30 x 30 lattice under the 
periodic boundary condition. We explicitly checked for three sets of parameters 
{t\\/U, fx/U) marked as points of 1, 2 and 3 in Fig.[6]B. The stripe-ordered ground 
state with a 2 x 2 unit cell depicted in Fig. |4] A is found stable against small 
random perturbations in all the three cases. Then we further apply the Gutzwiller 
type theory assuming a 2 x 2 unit cell and obtain the phase diagram depicted in Fig. 
[S]A which includes both the stripe ordered superfluid phase and Mott-insulating 
phases. 

In order to gain a better understanding of the numerical results, we write the 
trial condensate with the p-orbital configuration on each site as 

e'^'^^{cosa\px) + iaf!S\na\py)). (28) 

We have checked that the optimal pattern for the U{\) phase 0^ does not depend 
on the orbital mixing angle of a, and it also remains the same for all the coupling 
strength. The phase mismatch defined in Eq.[27|on the tilted bonds for a general 
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orbital mixing angle a can be calculated through simple algebra as 

A6I = 27 - 7r/2, with tan 7 = Vii&ria, (29) 

and the corresponding Josephson current is j — nt sin A9. The value of a is de- 
termined by the minimization of the energy per particle of the trial condensate 
as 

TT TlU TliU 

E(a) ^ -i[l + 2sin(2a+ -)] - —sin^2a + —, (30) 

6 6 3 

where the first term is the contribution from the kinetic energy which requires phase 

coherence, and the second term is the interaction contribution reflecting the Hund's 

rule physics. In the strong and weak coupling limits, the energy minimum is located 

at a = ^ and ^, respectively. The corresponding fluxes in each rhombic plaquette 

(f> = 4A6'/(27r) — and ±i respectively, which agree with the previous analyses. In 

the intermediate coupling regime, we present both results of $ at rt = 3 based on 

the Gutzwiller mean field theory and those of Eq. ([50]) in Fig. [SB. They agree with 

each other very well, and confirm the validity of the trial condensate wavefunction. 

Moreover, in the momentum space, the trial condensate for a general a can be 

expressed as | -i= (■0^^ + i^K:,)\ '^here V'ks.sI^) = e'-^''''*'l02,3(a)) with 

\4>2,3{oi-)) — —cosa\px) ^suia\py) respectively. 

2.6. Orbital angular momentum ordering in the Matt-insulating 
states 

So far, we have only discussed the orbital ordering in the superfluid states. Since 
the 0AM moment is a different degree of freedom from the superfluid phase, we 
expect that its ordering can survive even inside the Mott-insulating phases. In this 
subsection, we continue to study the exchange physics of orbital bosons and the 
related ordering of 0AM moments in the absence of the superfluidity order. For 
simplicity, we will use the triangular lattice as an example. 

We consider the Mott-insulating phases with n spinless bosons per-site and two 
degenerate orbitals of and Py. We define the TR doublets of all the particles in 
the states of px ± ipy as the eigenstates of the Ising operator az with the eigenvalues 
of ±1. The Ising part of the effective exchange Hamilton occurs at the level of the 
second order perturbation theory, while the Ising variable fiipping process occurs 
at the 2n-th order perturbation theory. We consider the large-n case in which the 
system is deeply inside the Ising anisotropy class, and no orbital-fiip process occurs 
at the leading order. In the following, we will study the physics of the 2,3,4-site 
exchange processes. 

2.6.1. The two-site exchange 

Let us consider the virtual hopping processes in the Mott-insulating states along 
the bond depicted in Fig. [5] A. Both the cr-bonding and 7r-bonding t± are kept. 
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With the definition of pj_ = -^(j>\ i ipp, we can express the hopping as 

aa' 

- ^ E + h.c.} + Y: {pI.P.-^^-''^'^ + h.c.}, (31) 



where the definition of angles 9i, 9j foUows the convention depicted in Fig. [5] A. 

We calculate the energy shifts in both the ferro-orbital configurations (di = 
(Tj) and the antiferro-orbital configurations { Oi = aj) within the second order 
perturbation theory. The energy difference between these two configurations is 

2n(n+l)[(t|| -tx)/2]^ 2n[(t|| + t^)/2]^ 
AE[ AE2 

2n(n+l)[(t|| +tx)/2]^ 2n[ft|| - tj^)/2]^ 
= AEl AE, ' (^^^ 

where AEi — and AE2 = |C/(ri + 1). The antiferro-orbital configuration has 
lower energy, which arises from the fact that the orbital-flip process has a larger 
amplitude than that of the orbital non-flip process described in Eq. [211 This is 
because that the cr-bonding and vr-bonding amplitudes have a 7r-phase shift. The 
Ising coupling Jafo reads 

H = Jafo E '^S)^z{3) (33) 



where Jafo ~ ^"3;/ *" with only the leading order contribution kept. If t± is set 
to zero, Jafo varnishes. This is clear from the fact that we can flip the sign of the 
p-orbit component perpendicular to the bond direction. This operation changes the 
value of az, but has no effect on the energy. In this case, we need to further study 
the multi-site virtual hopping processes. 



2.6.2. The three-site ring exchange 

The (T-bonding term by itself gives non-zero ring exchange terms for the multi- 
site processes, thus we neglect the contribution from the 7r-bonding part. In the 
following, we only keep the leading order virtual process proportional to (nt||)^/C7^. 

We consider a triangular plaquette with three sites {i,j, k) each of which is de- 
noted by the particle number and the Ising variables as (n, Uz^i)-, (ji, (Jzj), {n, <Jz,k), 
respectively. There are 12 different ring-hopping processes whose contribution is at 
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the order of {nt^Y' /U'^ . We enumerate 4 of them exphcitly as 

(n,cr^,i)(n,cr2j)(n, cr^^fe) {n - l,(T^,i)(n + l,az,]){n,(T^^k) 
^ {n- l,o-2,i)(n, cr2j)(n + l,(T^,fe) {n,(Tz,i)in,(Tz,j){n,az^k) 
{n,az.i)in,(Tz.j){n,az^k) ^ {n + l,az,i){n- 1, (T2j)(n, cr^^fc) 
{n,az.i)in- l,cr2j)(n + l,(Tz,fe) {n,az,i)in,(Tz,j){n,az^k) 
{n,az,i)in,az,j){n,crz^k) {n - l,cr^,i)(n + l,az,]){n,(Tz,k) 
{n,(7z,i){n+ l,o-^j)(n- l,cr^,fe) ^ {n,az^i){n,azj){n,azM) 
{n,az,i)in,(Tz,j){n,az,k) (n + 1, (T2,i)(n - l,az^j){n,az,k) 
^ {n~ l,az,i){n,az,j)in + l,(T^,fc) ^ {n,(Jz.i){n,(7z.j){n,(Jz,k)- (34) 

The other 8 processes can be obtained by a cyche permutation i j k. After 
sum over all these processes, the corresponding energy shift reads 

Ai?(a,,a, ,a,)«6(^fKn + l))3/2f:!^^^t^ = - J3 cos(2^$3), (35) 

where J3 — 2{^*Ey ^"^^ = $3 is defined as before in the superfluid case. 
Again for each plaquette, $3 takes the value of ±i which correspond to azS taking 
the values of two Is and one —1 or two —Is and one 1. Eq. 1351 plavs a similar role 
to vortex core energy in Eq.[2H Eq.[35]by itself has a huge ground state degeneracy 
because the requirement is the same as that of the antiferromagnetic Ising model. 
We need to add the four-site process to lift the degeneracy. 



2.6.3. The four-site ring exchange 

The four-site process mimics the interaction between two adjacent vortices. By 
the spirit of perturbation theory, the 4-particle process is further suppressed by 
a factor thus the interaction between vortices should be of short range. This 
makes sense in the Mott insulating state due to the condensation of vortices and 
the resulting screening effect. 

Again we only consider the effect from the cr-bonding term, and only keep the 
leading order contribution at the level of {nt)'^/U^. Similarly, we have 48 different 
virtual hopping processes. For a plaquette with vertices {ijkl), the energy shift is 

AE{ai,aj,ak,(Ji) « - J4 cos(27r$4), (36) 

where J4 = 3(^ae)^ ■ ^^"^ square lattice, $4 = (ai + (T2 + (73 + 0-4) /4. If ij^ = 0, Eq. 
[Ml is the leading order contribution, thus the effective Hamiltonian looks the same 
as that in the superfluid case but with a reduced coupling constant. If ij^ 7^ 0, then 
the two-site exchange of Eq. [33] results in the long range antiferro-orbital order. 

In the triangular lattice, the <I>4 of the four-site plaquette (ijkl) just equals the 
sum of the flux of the two triangular plaquettes. Since each triangular plaquette 
takes the flux value of ±i, the rhombic four-site plaquette can only take the flux 
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Fig. 7. A)The radial density distribution of spin up and down components, and the total density 
distribution in the unit of A^o at a = 4 and /3 = 40. B) The spin density distribution along the 
rr-axis which spirals in the z-x plane at an approximate wavevector of 2ko, whose value at the 
origin is normalized to 1. The spin density in the whole plane exhibits the skyrmion configuration. 
From Wu et al. Ref. 5. 



value of 0, ±i. The stripe ordering pattern gives the maximum number of the 
zero-flux four-site plaquettes, and thus is selected as the ground state. The stripe 
ordering appears at temperatures T < J4 and disappears at J3 > T > J4 in which 
the three-site exchange process dominates the physics. 

3. Exotic condensation of bosons with spin-orbit coupling 

In this section, we discuss another route for "complex-condensation" of bosons, i.e., 
spinful bosons with spin-orbit (SO) coupling. The proof of Feynman's "no-node" 
theorem crucially relies on the fact that the kinetic energy of bosons only contains 
even power of momentum, thus SO coupling which linearly depends on momentum 
invalidates Feynman's proof. In the following, we consider the simplest example of 
the two-component bosons with the Rashba-like SO coupling and investigate their 
condensate in the inhomogeneous harmonic traps. A related work on the BEG of 
two-component bosons has also been studied by Stanescu et al. in Ref. 

We consider bosons with two-components denoted as pseudospin up and down. 
The many-body Hamiltonian for interacting bosons is written as 

+ hXaipli-iS/ ya^ + i\/x(yy)i'() + ^'4'ii^l'4'0ipa, (37) 

where tpa is the boson operator; a refers to boson pseudospin t and J,; Vext is the 
external potential; g describes the s-wave scattering interaction which is assumed 
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to be spin-independent for simplicity; Ai?, is the Rashba SO coupling strength. Al- 
though Eq. [37]is of bosons, it satisfies a suitably defined TR symmetry like fermion 
systems as T = ia2C with = — 1 where C is the complex conjugate and 02 
operates on the pseudospin doublet. In the homogeneous system, the single parti- 
cle states are the helicity eigenstates of ct • (fc x z) with the dispersion relations of 
£±(^) = 517 (^"F^o)^ where fcg — ■ The energy minima are located at the lower 
branch along a ring with radius fco. The corresponding two-component wavefunc- 
tion '4>+(k) with |fc| = fco can be solved as V'+(^) = :^(e~*'^'=/^, ze*'^'=/^), where 4>k 
is the azimuth angle of k. 

Instead of discussing the condensation in the homogeneous system in which 
frustrations occur due to the degeneracy of the single particle ground states, we 
consider a 2D system with an external harmonic trap with Vexir) ~ ^Muj^r^ . We 
define the characteristic SO energy scale for the harmonic trap as Eso = fiXn/l 
where I = yjh/ [Mlot) is the length scale of the trap, and correspondingly a dimen- 
sionless parameter as a = Eso/ {^u^Jt) = Iko. The characteristic interaction energy 
scale is defined as Eint — gNo/{TTp) where Ng is the total particle number and the 
dimensionless parameter /? — Eint / {ft^Jr) ■ 

The Gross-Pitaevskii equation can be obtained as the saddle-point equation of 
Eq. [32] as 

= SVa(^,0)- (38) 



Due to the 2D rotational symmetry, the ground state condensate wavefunctions can 

1 ,1 

2 2 ^ 2' 



be denoted by the total angular momentum = + ^(Tz = which can be 



represented in the polar coordinate as 



where r and are the radius and the azimuthal angle, respectively. Both /(r) 
and g{r) are real functions. Eq. [35] is numerically solved in the harmonic trap 
with parameters of a = 4 and (3 — 40. We choose the condensate as one of the 
TR doublet a-nd present its radial density profiles of both spin components 

|/(r)p and |5(r)p at a = 4 in Fig. [7] A. Further in this strong SO coupling case 
with a ^ 1, the condensate wavefunction has nearly equal weight in the spin up 
and down components, i.e. J drd(j) r\f{r)\'^ « J drdcj) r|g(r)p, thus the average spin 
moment along the z-axis equals to zero. The total angular momentum per particle 
jz = f is mainly from the orbital angular momentum polarization, i.e., one spin 
component stays in the s-state and the other one in th e Px + ipy-state. This is an 
example of half-quantum vortex configuration ISZESIM]^ th^s spontanously breaking 
TR symmetry. Clearly this is a "complex- valued" ground state wavefunction beyond 
the "no-node" theorem. 
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This condensate wavefunction exhibits interesting spin density distributions in 
real space as skyrmion-Hke spin textures. The radial wavefunction in pseudo-spin 
up and down components f(r) and g{r) exhibit oscillations with an approximate 
wavevector of fco, which originates from the ring structure of the low energy states 
in momentum space, thus are analogous to Friedel oscillations in fermion systems. 
The pseudo-spin up component is s-wave like, thus /(r) reaches the maximum at 
r = 0; while the down component is of the p-wave, thus g{r) = at r = 0. In other 
words, approximately there is a relative phase shift of ^ between the oscillations 
of /(r) and g{r). In Fig. [7] A, and are plotted. The spin density 

distribution 5(r, 0) — ip\ (r, (j>)aaaipi a{r, 4>) can be expanded as 

S.{r.4>) = \{\f{r)? - IffMH, 5,(r,0) = /(r)g(r) cos 0, 

Sy{r,4,)^ f{r)g(r)sinc^. (40) 

Along the x-axis, the spin density lies in the z-x plane as depicted in Fig. [7| B. 
Because the circulating supercurrent is along the tangential direction, the spin den- 
sity distribution is along the radial direction and exhibits an interesting topological 
texture configuration which spirals in the z-x plane at the pitch value of the density 
oscillations. The distribution in the whole space can be obtained through a rotation 
around the z-axis. This spin texture configuration is of the skyrmion-like. 

Next we discuss the possible realization of the above exotic BEC in exciton 
syste ms. Excitons are coin posite objects between conduction electrons and valence 
holes Iini 41 | 42 | 43 | 44 | 45 | 46 | jj^ particular, t he recently p rogress on the indirect exci- 
ton systems greatly enhances the lift-time I 47|48|49 |50j^ which provides a wonderful 
opp ortunity to investigate the exotic state of matter of the exciton condensation 

The ordinary bosons are too heavy to exhibit the relativistic SO coupling in 
their center-of-motion. Due to the small effective mass of excitons, SO coupli ng can 
result in important consequences, including anisotropic electron-hole pairing ^^2221^ 
spin Hall effect of the center-of-mass motion of excitons ^"^ and the Berry phase 
effect on exciton condensation 

We will show that the Rashba SO coupling in the electron band also survives in 
the center-of-mass motion of excitions. We begin with the Hamiltonian of indirect 
excitons 

^2 

He = + dly) + ih\R^el{de,x<yy " de,y<7x), 



2m^ 
2m 



hh 

^"^^^ " /l^ - 12 ^^-^^ 

where m* is the effective mass of conduction electrons; X^^ei is the Rashba SO 
coupling strength of the conduction electron, e is the dielectric constant; d is the 
thickness of the barrier. 
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For small exciton concentrations, we only need to consider the heavy hole [hh) 
band with the effective mass and jz = ±|?i, which is separated from the light 
hole band with a gap of the order of 10 meV. We consider the center-of-mass motion 
in the BEG limit of excitons, wh ich can be separated from the relative motion in 
He and Huh- Similarly to Ref. the effective Hamiltonian of the 4-component hh 
excitons denoted as {se,jhh) — (i^j f )i (ij; ) can be represented by the matrix 
form as 

I Eex{k) Hso{k) \ 

^ ^ H*so{k) Eex{k) + A{k) W{k) 

W*(k) Eexik) + A{k) Hsoik) ' 

V HUk) Eexik)/ 

where k is the center-of-mass momentum; Af = m* + m'^i^ is the total mass of the 
exciton; Eex{k) = h?k'^ /{2M); Hso{k) = —^XR,ei{ky +ikx); A(fc) is the exchange 
integral and W{k) has the d-wave structure as {kx + iky)^, both of which are 
exponentially suppressed by the tunneling barrier for indirect excitons, and will 
be neglected below. Consequentially, Hex becomes block-diagonalized. We consider 
using circularly polarized light to pump the exciton of (— ^, f ), and then focus on 
the left-up 2x2 block of the Hex matrix with heavy hole spin jz — ^- If we use the 
electron spin number "f and J, as the exciton component, we will arrive the Eq. [37] 
with the renormalized SO coupling strength Xn = \ji^ei'm*e/M . 

We next justify the above choice of the values of parameters based on exper- 
imental situations. The effective masses of electrons and holes i n GaAs/AlGaAs 
quantum wells are m* ~ 0.07me and mj^ w O.lSme in Ref. I ^^ l ^^ l -pj^g fj^ashba SO 
coupling strength hXa can reach 1.8 x 10"^^ eV.m in Ref.^^. Thus we can estimate 
a reasonable value of fco = "~^r^ ~ 1-6 x 10** cm~^. For a harmonic trap with 
I = 2.5/xm, a — k^l w 4 and hior = h'^ /(MP) = 0.5 niK. In two dimensional har- 
monic traps, the critical condensation temperature Tc ~ hcoTy/No^^- If we take the 
exciton density p = 5 x 10^° cm~^ and the effective area tt P, w e arrive at « 50 
mK, which is an experimentally available temperature scale The average inter- 
action energy per exc iton in the typical density regime of 10^*^ cm^^ is estimated 
around 2meV in Ref. thus we take the interaction parameter /3 = 40 in the 
calculation above. The spatial periodicity of the spin texture is about 7r/fco ~ 2fj,m 
and, thus, is detectable by using optical methods. 

4. Conclusion 

We have reviewed the exotic condensations of bosons whose many-body wavefunc- 
tions are complex-valued in the coordinate representation, thus they go beyond 
the well-known paradigm of the "no-node" theorem. We studied two possible ways 
to bypass this theorem to achieve unconventional condensations, i.e., meta-stable 
states of bosons in the high orbital bands in optical lattices, and spinful bosons 
with SO coupling. 
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The first mechanism of orbital bosons is essentially an interaction effect which 
is characterized by the Himd's rule. With the orbital degeneracy, bosons favor to 
enlarge their spatial extension to reduce the inter-particle repulsion, which results 
in the maximization of their onsite 0AM. We reviewed the ordering of the 0AM 
moments in both the square and triangular lattices due to the inter-site coupling 
in both the weak and strong coupling limits. The low energy excitations include 
both the gapless phonon mode and the gapped orbital-flip mode. The survival of 
the 0AM ordering in the soft Mott-insulating regime is also discussed. The second 
mechanism employing SO coupling, which linearly depends on momentum, is a 
kinetic energy effect. In this case, taking the absolute value of a non positive- 
definite wavefunction will change its energy and thus invalidate Feynman's proof. 
We have shown that the condensate wavefunction in a harmonic trap becomes a 
half-quantum vortex and develops skyrmion-like spin textures. In both cases, TR 
symmetry is spontaneously broken which are not possible in the conventional BEG. 

There are a number of interesting open issues for future study. For example, 
the ordering of 0AM moments in various different lattices and in three d ime nsions 
generates a variety of challenging problems of the lattice gauge theory More 
importantly, in order to facilitate the communication between the theory work 
with cold atom experiments, detailed calculation and optimization of the life time 
of orbital bosons in different lattices are highly desired. Furthermore, the dynamics 
of bosons in high orbital bands in another challenging topic. Taking into account 
the tremendous progress of cold atom physics, it would be great if this research 
can stimulate general interests on the unconventional condensates beyond the "no- 
node" theorem. 
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